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ABSTRACT 

This paper reports two possible corotating shocks in the inner helio- 
sphere where the solar wind is composed of low-6 plasma. In the region where 
the solar wind is slightly super-Alfvenic , reverse corotating fast MHD shocks 
can form at the leading edge of a high-speed stream. These shocks possess a 
switch-on mechanism for amplification of tangential small-scale fluctuations. 

The second one is the coronal slow shocks which may be imbedded in large coro- 
nal holes at low altitude in the sub-Alfvenic region. 

Introduction 

In a heliocentric coordinate system corotating with the sun, the large- 
time-scale structure of the solar wind is often quasi-steady for a few solar 
rotations. Many large-scale structures of the global solar wind may be treated 
as steady-state field-aligned flows^in a frame of reference corotating with the 
sun at a const an t angular velocity ft. The flow velocity in the corotating frame 
U is either parallel or anti-parallel to the magnetic field The solar wind 
process may be considered to consist of an expansion of the solar wind along 
each corotating stream tube and a dynamical interaction between neighboring 
stream tubes. 

Recently, magnetohydrodynamic (MHD) treatments have been employed to study 
the interaction of corotating high-speed streams in the highly super-Alfvenic 
region of the heliosphere by several authors. For example, Pizzo (1982) has 
studied some 3D problems under the adiabatic assumption, and Whang and Chien 
(1981) have studied some MHD interaction problems near the solar equatorial 
plane to include the heat conduction process. 

This paper reports two possible corotating MHD shocks in the low-B region 
of the inner heliosphere. The solar wind flow is sub-Alfvenic near the sun, 
and a transition from sub-Alfvenic to super-Alfvenic takes place at a heliocen- 
tric distance of approximately 0.1 AU. In the slightly super-Alfvenic region, we 
studied the interaction between two streams of solar wind at different Alfven 
numbers, a slow stream at a lower Alfven number precedes a fast stream. As a 
result of interaction, reverse corotating fast MHD shocks may form near 0.1 AU 
at the leading edge of the high-speed stream. These shocks possess a "switch- 
on" mechanism for amplification of tangential small-scale fluctuations. The 
second shock is the coronal slow shock, which may be imbedded in large coronal 
holes at altitudes of a few solar radii in the sub-Alfvenic region. 
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Mathematical Development 


The governing equations for the expansion of the solar wind along each 
stream tube may be represented by a system of five equations. The equation 
of continuity is 


e • V (pU/B) = 0 

s 


( 1 ) 


where p is the mass density, and e = i?/U is the unit vector along the stream- 
line direction. The parallel component of the equation of motion is 



(V — j + ~ V p - H) 


0 


( 2 ) 


where p is the thermal pressure, and 


H = -(GM/r 3 ) r — 2 ft x U - Q x (3 x r) 


G is the gravitational constant, M is mass of the sun, and r the heliocentric 
position vector. The force acting on a unit mass of the solar wind plasma, S, 
consists of the gravitational force and two fictitious body forces: the 
coriolis force and the centrifugal force. 


The variation of p and p along each streamline is governed by 

e g • (Vp - c 2 V p) = 0 (3) 

h 

where c = (ap/p) is the sound speed in a heat conducting plasma, a is the 
polytropic index. The conduction heat flux is parallel to e g , 

q = -IC ®s • v T (4) 


where T is the one-fluid temperature, and the thermal conductivity 
K - «I« 


with K - 7.5 x 10 erg s cm deg * . The polytropic index is directly 
related to the variation in thermal state and in heat conduction process along 
each streamline, 

n • { <! - “> 7 » - lio 7 <f> } - 0 <5 > 

where m is the proton mass and k Boltzmann’s constant, a is not a constant along 
each streamline. In the expansion of the solar wind, a < 5/3 when conduction 
heat acts as an energy source, it continuously converts into other forms of 
energy, ot = 5/3 in the region where q is a divergence-free vector (including 
the adiabatic condition with q = 0) . A special feature of the present formula- 
tion of energetics for a heat conducting plasma is the use of a varying poly- 
tropic index a in (3) and (5) instead of the equivalent equation of energy con- 
servation 
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From equations (1) - (3), we can obtain 
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0 


( 6 ) 

(7) 

( 8 ) 


-e • V £n B measures the rate of fractional change of cross-sectional area of 
a stream tube along the streamline direction. Once e s and the variation of 
cross-sectional area are known, the expansion equations can be integrated along 
each streamline. 


The variation in the direction and the magnitude of the cross-sectional 
area of each stream tube is controlled by the interaction between neighboring 
stream tubes (Whang, 1980). The interaction process is governed by the 
divergence-free condition of the magnetic field. 

V • B - 0 (9) 

and the dynamical equilibrium of momentum transverse to the streamline 
direction 


(U 2 - a 2 ) e • V e + (f - e e ) • C— V d* - H) = 0 f 10) 

s s s s p ' 

where a = B/(4irp) 2 is the Alfven speed, I is a unit dyadic, and p* = p+B 2 /8tt 
is the sum of thermal and magnetic pressures. Making use of equation (6), we 
can write (9) in a new form 

£u 2 (c 2 +a 2 ) - c 2 a 2 j V • e g +e s * ^ (U 2 -c 2 ) ^ V p*+c 2 H J = 0 (11) 

The system of the interaction equations, (10) and (11), shows that the varia- 
tion of the streamline direction is controlled by Vp* and H. The system is 
elliptic or hyperbolic, depending on the sign of 

[u 2 (a 2 + c 2 ) - a 2 c 2 ] / (U 2 - a 2 )(U 2 - c 2 ). It is hyperbolic in the super- 
Alfvenic, supersonic regime 

U > max (a, c) 


or in a sub-Alf venic , subsonic regime with 

( — + — ) < U < min (a, c) 

V a 2 c 2 / 


The system of interaction equations is elliptic in other regimes. 
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In the two hyperbolic flow regimes, solutions of the interaction equations 
respectively represent families of fast and slow MHD waves. When successive 
MHD waves travel faster than the preceding ones, they overtake and combine with 
each other. This phenomenon leads to the formation of fast or slow MHD shocks. 

Corotating Shocks Near 0.1 AU 

A simultaneous solution of the expansion equations and the interaction 
equations is carried out to demonstrate the formation of corotating shock near 
the equatorial plane in the hyperbolic flow region outside -0.1 AU where the 
flow of a low-3 plasma is supersonic and slightly super-Alfvenic. We assume 
that near the solar equatorial plane the flow properties are independent of the 

latitude. Let <f> denote the directional angle between -r and e , e = 5^ x e /Q 

^ s’ n s 

a unit vector normal to the streamline direction, H = e • H and H = e • S. 

s s n n 

Then we may write the interaction equations, (10) and (11), near the equatorial 
plane as 

(U 2 - a 2 ) ( e • V <j> - ■ —■ ■ * ) + - e • V p* = H (12) 

s r p n n 

and 

[u 2 (c 2 + a 2 ) - c 2 a 2 ] (I n • V * + ) 

+ i (U 2 - c 2 ) t • 7 p* - -c 2 H (13) 

p s s 

In the super-Alfvenic hyperbolic regime, the Mach angle y for fast MHD 
waves is defined by 

cos y = (1 - a 2 /U 2 )^ (1 - cVu 2 )* 5 (14) 

We may introduce 

-)■ 

c, = cos y e + sin y e , (15) 

± s n 

unit vector^ along the two characteristic directions. The characteristic curve 
defined by c + running ahead of the streamline in the direction of solar rota- 
tion may be called the forward characteristic curve, and that defined by c_ 

running behind the streamline the reverse characteristic curve. We can 
express the interaction equations in the form of characteristic equations 
obtained from a linear combination of (12) and (13), 

c ± • (V <f> ± Q V p*) = S + (16) 

where Q = — 

P (U 2 - a 2 ) tan y 

cos y H ___ 

and = — + 

1 U 2 - a 2 


c 2 H 


U 4 sin y 


H 


sin (<t>±y) + sin 


■] 
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Solutions for MHD interaction between solar wind streams in the low-g 
region of the inner heliosphere can be obtained from integrations of the expan- 
sion equations along streamlines, and integrations of the interaction equations 
(16) along characteristic curves. A method of solution for solar wind streams 
in the super-Alfvenic region has been developed previously (Whang and Chien, 
1981). The present approach differs from the previous one in two aspects: 

(i) using a varying polytropic index in (3) and (5) instead of the equivalent 
equation of energy conservation, and (ii) introducing a new form (11) to re- 
present the divergence-free condition (9) . The advantages of the present 
approach are (i) each of the final interaction equations (16) contains only 
derivatives along a characteristic direction, and (ii) slow mode MHD waves are 
not suppressed. 

The following numerical solution demonstrates the interaction between two 
streams of solar wind flow at different Alfven numbers, in which a slow stream 
precedes a sharply bounded fast stream of higher Alfven number. Assume that 
at the initial surface r = 22.5 R , the profile for the radial speed (in km/s) 
be given by 


U = 400 + 8 (1 - erf w) 

r 

where w is the azimuthal angle in degrees. A smooth increase of the flow speed 

from U = 400 km/s and U/a = 1.087 to U = 416 km/s and U/a = 1.132 takes place 
r r 

in a narrow region with its longitudinal width of approximately 4°. We further 
assume that at the initial surface 

U = - fir 

0) 

“2 — 1 

nU^ = 2,795 x 10 5 protons cm s 
and T = 4.3 x IQ 5 K 

From the slow to the fast stream, B changes from 484 to 439 nT, and p* from 
1.016 x 10" 6 to 0.847 x 10 _6 dyne/cm 2 . 

We have obtained the numerical solution for all flow parameters in the 
region 22.5 £ r/R < 30.4. In the top frame of Figure 1 are plotted the 

streamlines (in heavy curves) and the two families of characteristic curves in 
the physical plane. The streamline passing through w = 0 at the initial sur- 
face is identified as the central streamline. At the initial surface, the Mach 
angle decreases across the interface from y = 83.6° in the slow stream to 
y = 62.8° in the fast stream due to changes in Alfven numbers. Thus the slopes 
of those characteristic curves starting near the initial interface change rap- 
idly. The merging of the reverse characteristic curves shown in the plot 
indicates the formation of a reverse corotating shock, it signifies a piling up 
of the fast-mode MHD waves to form a fast MHD shock. The lower frame of Figure 
1 shows the profile of the total pressure p* at various heliocentric distances. 
Shock fronts with steep profiles of p* appear at the leading edge region. The 
shock fronts become broader and weaker at approximately 15° (or ~24 hours) 
behind the central streamline due to geometrical effects. Unlike the total 
pressure, other flow properties (U, p, T, etc.) change abruptly across the 
stream interface around the central streamline. Disturbances of flow proper- 
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perties due to stream interaction are confined to an interaction region bounded 
by the two characteristics which intersect the initial surface near the central 
streamline. The solar wind streams undergo a simple expansion process without 
strong interaction between neighboring stream tubes outside of the interaction 
region. In the region where the flow is slightly super-Alfvenic, one can distin- 
guish clearly the shock front from the interface because of large Mach angles. 
Between the interface and the reverse shock, the flow properties carry the 
signatures of a shocked plasma as shown in Figure 2, which plots the profiles 
of U r and T at various heliocentric distances. Although the temperature is 

assumed to be uniform at the initial surface, the shocked plasma has a higher 
temperature in the leading edge region. 

This example illustrates that in the low-g region of the inner heliosphere 
where the solar wind flow is slightly super-Alfvenic, if the Alfven number in- 
creases at the leading edge region of a high-speed stream even at a very slow 
rate, reverse corotating shocks may form in the region. These shocks are fast 
MHD shocks of low-g plasma, with the upstream normal Mach number slightly great- 
er than 1. They are capable of amplifying the tangential component of the mag- 
netic field. Thus, the presence of a small random fluctuation field upstream 
of these shock surfaces can be amplified to produce a large, tangential fluc- 
tuation field downstream of the shock. This effect may be called a switch-on 
mechanism, which will be discussed in detail in a separate paper. This ampli- 
fication mechanism is probably responsible for the presence of large-amplitude 
anisotropic fluctuations observed at the leading edge of high-speed streams 
(Belcher and Davis, 1971). 


Coronal Slow Shocks 

At altitudes of a few solar radii, the solar wind is a sub-Alfvenic low-g 
plasma. It has been suggested that, in this region of the inner heliosphere, 
corotating slow MHD shocks might be imbedded within large polar coronal holes 

/m i noo\ 

V Writing , x yoz. ) • 

Figure 3 depicts the geometry of a coronal slow shock around a simple two- 
hole corona. The solar wind originating from open-field coronal holes is ac- 
celerated to high speeds at low altitudes in the corona. Streams originating 
from the edge of the polar open-field regions flow around the curved boundary 
of the helmet-shaped, closed-field region (Pneuman and Kopp, 1971). Near the 
inner edge of the neutral sheet, the flow direction of each stream changes 
suddenly, becoming parallel to the neutral sheet. An oblique slow shock can 
develop near or at the neutral point. The shock extends polewards to form a 
composite surface of discontinuity surrounding the sun. These standing shocks 
could occur quite close to the solar surface, at 2 to 5 solar radii above large 
coronal holes. Numerical solutions have been constructed for the expansion of 
the solar wind in open field region of the solar corona along stream tubes with 
prescribed geometry by Pneuman and Kopp (1971) and by Whang (1983). These 
solutions did not consider dynamical interaction between neighboring stream 
tubes, and their streamline geometry was introduced under an a priori assump- 
tion that no shocks exist in the corona. 

For not very large coronal holes, the coronal slow shocks are weaker at 
lower altitudes, and could have a partially concave surface as shown in Figure 4. 
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figure 4. For not very large coronal holes, the slow shocks are weaker at lower 
altitudes, and could have a partially concave shock surface. 
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Let 0 denote the angle between the magnetic field and the shock normal, 
and C g be the normal slow speed, 

C 2 = (a 2 + c 2 ) - | ^(a 2 + c 2 ) - 4a 2 c 2 cos^G ^ 

The normal component of the flow speed U is greater than C upstream, and less 

n s 

than C downstream of a slow MHD shock. The flow downstream of a slow shock is 
s 

hyperbolic, sub-Alfvenic and subsonic. In the hyperbolic flow region, solu- 
tions of the interaction equations (10) and (11) represent (obtuse) Mach angle 
y with 


cos p 



(17) 


(Sears and Resler, 1961). Successive slow waves travel upstream faster than 
the preceding ones, so that they overtake and combine with each other to form 
a slow shock upstream of the sub-Alfvenic and subsonic hyperbolic flow region. 


In contrast with the jump conditions for fast MHD shocks, B and 0 decrease 
across a slow shock. $ increases due to the sudden rise in the thermal pressure 
and the drop in the magnetic pressure. However, p*, as well as p and T, always 
increase. The possible existence of coronal slow shocks can significantly 
revise our understanding on the dynamical structure of the global solar wind. 
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